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In their _letter on the quantum Child-Langmuir law, 
Ang et al [l] include exchange correlation effects within 
the Kohn-Sham density functional theory (DFT) [2] and 
explore numerically the maximum transmitted current 
{J max) in nanogaps. We show here that the calculations 
are in error as the exchange-correlation component of 
the chemical potential has been ignored while fixing the 
boundary conditions for the Hartree potential. 

The analysis in [IJ is based on solving the time- 
independnent Schrodinger equation —(P^p/dx^ + Vef fifj — 
Etp with an effective potential energy, Ve// = —eV + 
Vxc X Eh, where Eh = e^/(47reoao) is the Hartree en- 
ergy, ao the Bohr radius, V^c = ^xc - {rs/^)dea:c/drs, 
Ts = [3/(4?™)]-'^/'^ is the Wigner-Seitz radius, n the 
electron number density and e the magnitude of the 
electronic charge. In the above, e^c is the exchange- 
correlation energy density within the local density ap- 
proximation, and V is the Hartree potential satisfying 
the Poisson equation cPV/dx^ = en{x)/eo = e|?/'(a;)P/eo 
with boundary conditions V{Q) = and V{D) — Vg 
where Vg is the gap voltage difference and D the gap size. 
These equations are solved in a nanogap with the wave- 
function and its derivative matched at the collector 
boundary under certain assumptions (Ref. [sj addresses 
one of these). 

We first note that the applied voltage difference, 
Vg — —{nc — ^'e)/g where and ^jle refer respec- 
tively to the chemical potential at the collector and in- 
jection planes. For convenience, we consider the ref- 
erence as — — eV^(O) -I- T4c(0) x Eh = so that 
E = Q refers to injection from the Fermi level. Thus 
V(Gi) — Vxc{^) X En/e. It follows that the chemical po- 
tential at the collector is ~eV{D) + Vxc{D)xEH = -eVg. 
Thus V{D) ^Vg + Vxc{D) xEH/e. In writing the above, 
we have implicitly assumed continuity of the chemical 
potential at the interfaces under steady-state conditions. 
Note that when exchange-correlation is neglected alto- 
gether, the boundary conditions for V are 1^(0) = and 
V{D) = Vg respectively as assumed in [ij. 

When Vxc is substantial, and the boundary conditions 
chosen are V{D) = Vg and V^(0) = 0, the results can be 
unphysical. As an example, consider the case D — Inm, 
Vn = O.IV and E = Q. Using the formalism of Ref. 
Ul, Jmax turns out to be Jmax — 278 JcL where Jcl is 
the classical Child-Langmuir current density. The cor- 
responding effective potential energy is shown in Fig. [T] 
Note that Fe// is negative everywhere. The injection en- 
ergy E has clearly no relation to the effective potential 
at either end of the nanogap. 

The corrected boundary conditions for V lead to 



FIG. 1: The effective potential energy Veff — —cV + Eh x V^c 
for D = Inm, Vg = O.IV (top) V"(0) = 0, V{D) = Vg and 
J^a^ ~ 27SJcL (bottom) V{0) = Kc(0) x E„/e, V{D) = 
Vg + Vxc{D) X Eh/c and Jmax ^ SAJcL- Tlie boundary 
conditions for tp axe the same as in Ref. 1] in both cases. 


Jmax — 3 A Jcl- The corresponding effective potential is 
shown in Fig. [T] (bottom). The error in determining the 
quantum Child-Langmuir law is therefore substantial. 

In regimes where Vxc is negligible (large Vg or D), 
V{D) — Vg and V{0) ~ are approximately the correct 
boundary conditions within the formalism of Ref. [l[ as 
we have verified for D = 50nm, Vg = 50V and E = 0. 

Finally, we note that the error in boundary conditions 
for V persists in subsequent publications by the authors 
of Ref. 1] . Results incorporating the rectified boundary 
conditions for V and improved boundary conditions for 
-0 will be published in a separate communication. 
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